Encoding many qubits in a rotor 
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We propose a scheme for encoding many qubits in a single rotor, that is, a continuous and 
periodic degree of freedom. A key feature of this scheme is its ability to manipulate and entangle 
the encoded qubits with a single operation on the system. We also show, using quantum error- 
correcting codes, how to protect the qubits against small errors in angular position and momentum 
which may affect the rotor. We then discuss the feasibility of this scheme and suggest several 
candidates for its implementation. The proposed scheme is immediately generalizable to qudits of 
any finite dimension. 



I. INTRODUCTION 

Quantum information processing, quantum simula- 
tion, and quantum communication require the individual 
manipulation and the coupling of quantum bits (qubits) . 
Early candidates for qubits were two-dimensional physi- 
cal systems, like the polarization of a photon or the spin 
of an electron, in which the state of the system directly 
represents the state of the qubit. Later, higher dimen- 
sional systems like atoms have been considered. There, a 
qubit is encoded in two energy levels while the remaining 
dimensions are ignored Infinite dimensional systems 
have been studied too. For instance, the encoding of a 
single qubit (actually qudit, a d-dimensional version of 
a qubit) in continuous variables and harmonic oscillators 
has been described in ^ and [s*] . Importantly, a common 
feature of all these encodings is the use of a single degree 
of freedom to encode a single qubit. 

An interesting alternative is to use a single degree of 
freedom to encode many qubits. Let us consider an infi- 
nite dimensional system which can be decomposed into a 
qubit subsystem and a system isomorphic to the original 
system itself, a feature only possible for infinite dimen- 
sional systems. Then the same decomposition can be 
repeated over and over again to obtain as many qubits 
as desired. This is the core idea of this article. We show 
how to encode N qubits (qudits) in a rotor, that is, a 
single continuous and periodic degree of freedom. The N 
qubits can be individually manipulated and coupled. Re- 
markably, the advantage of this approach is that a single 
unitary transformation acting on the rotor can simulta- 
neously rotate and entangle many qubits. Limitations 
will, however, inevitably appear with the physical imple- 
mentation. Furthermore, an encoding is practicable only 
if it can tolerate small errors affecting the physical sys- 
tem. Therefore, we use the stabilizer formalism [3, ^ to 
construct qubits that are robust against small errors in 
angular position and momentum. We then investigate 
practical approximations of our scheme and finally sug- 
gest a few implementations. 

In Sec.|TT]we briefly describe a quantum mechanical ro- 
tor and a qubit. We present in Sec. lIIIl a simple encoding 



to embed many qubits in a single rotor. We also give the 
corresponding universal set of quantum gates. In Sec. IIVI 
we provide a more general quantum error-correcting code 
to protect many qubits against small errors in angular 
position and momentum. It turns out that the protected 
qubits require an unlimited amount of energy. Thus we 
consider realistic approximations of the unphysical en- 
coded qubits in Sec.|Vl We then suggest in Sec. |Vl] pos- 
sible implementations in quantum optics, atom optics, 
and molecular physics. We briefly discuss the qudit case 
in Sec. EIIl We finally conclude in Sec. IVlIIl 



II. DESCRIPTIONS OF A ROTOR AND A 
QUBIT 

In quantum mechanics the angular momentum L for 
rotation around a fixed axis is a self-adjoint operator with 
a discrete and infinite spectrum. As we will see, this dis- 
creteness represents a very natural basis to encode qubits. 

We label the eigenstates of the operator L with an 
integer £ = 0, ±1, ±2, • • • as 



L\£)^£\£)^ 



(1) 



In this angular momentum basis {|^)}, the orthogonality 
and completeness relations take the simple form 



{£]£') = Se,e, and ^ \£){£\ = 1. 



(2) 



We further define the Fourier transform basis {16*}} of 
{\£}} as 



E 



(3) 



so that the two observables are complementary, that is, 

{9\£) = e'^^. (4) 

The angular position basis {\0)} is continuous and 
27r-periodic, \6 -t- 27r) = \6). Accordingly, we have 



\e') ^27Td'^^^\0-e') and / ^ 



1, (5) 
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where S^'^'^\(j)) is the 27r-periodic delta function defined 
^® '^kL-oo ~ 27rfc) and the integration is over any 27r 
interval. 

Let us now consider the operator e*"^, with real a. Its 
action on a ket \9) reads 

e»^|0) = |0-a). (6) 

The operator e*"^ is therefore called the shift operator 
in angular position. It follows from the fact that £ is 
an integer that e^"^ ~ 1. We proceed to introduce V, 
the shift operator in angular momentum. It is defined 
through its action on {|^)} as 

v\e)^\i + i). (7) 

Note for completeness that 

It is not too difficult to verify that these two shift oper- 
ators satisfy the commutation relation 

g^aL^ = T/e'"^e*". (9) 

These two properties together with e^*'^^ — 1 define the 
so-called Weyl pair or Schwinger operators of a rotor. 
The notion of a Weyl pair is closely related to the 
Heisenberg-Weyl and generalized Pauli groups [6,] . 

A quantum mechanical system is entirely characterized 
by its Weyl pair. First, each Weyl pair fully defines the 
Hilbert space of the system. Second, the Weyl pair is al- 
gebraically complete, that is, any operator acting on the 
system can be expressed in terms of these two operators 
only. This makes the Weyl pair extremely useful for dis- 
cussing the properties of a quantum mechanical system 
0i. 

As for the rotor we can introduce the Weyl pair {Z, X) 
of a qubit. The operators Z and X are both unitary and 
hermitian and such that 

Z^ ^X^ ^1 and ZX = -XZ. (10) 

By construction, Z and X are Fourier transforms of each 
other and each one corresponds to the shift operator on 
the eigenbasis of the other. 

Let us clarify here an important issue related to qubit 
encoding. For practical reasons, the involution property 
of the qubit Weyl pair might be withdrawn leading to 
the definition of a pseudo qubit. For example, one could 
only require Z^ and X^ to equal the projection operator 
onto the qubit subspace as seen in In this article, 
we focus our attention on genuine qubits for which the 
Hilbert space is the tensor product of the qubit space 
times the space for all other degrees of freedom. 

III. FROM A ROTOR TO MANY QUBITS 

For the sake of clarity we proceed in two steps. First, 
we show how to encode a single qubit on the basis of 
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TABLE 1: A qubit and a residual rotor. In this simplest 
case, the encoding corresponds to the parity of the angular 
momentum eigenstate The residual rotor is obtained by 
concatenating two adjacent angular momentum eigenstates. 

the angular momentum. This encoding can be schemat- 
ically summarized as rotor = qubit (8) rotor. Our map- 
ping corresponds to a splitting of the angular momen- 
tum space into two isomorphic subspaces. One sub- 
space corresponds to even angular momenta while the 
second subspace is associated with odd angular momenta. 
Such a decomposition is formally explained below but 
can be intuitively understood in terms of splitting the 
Fourier series of a 27r-periodic function into even and odd 
parts. The splitting effectively provides a qubit while 
each partial 27r-periodic Fourier series lead to a complete 
TT-periodic Fourier series, that is, a new rotor. Second, 
we use the same procedure and this remarkable tensor 
product structure to further split the remaining rotor. 
We can repeatedly factorize the remaining rotor degree 
of freedom and so create more and more qubit degrees of 
freedom. We shall now begin with a single qubit. 



A. Encoding a qubit in a rotor 

We want to encode a qubit subsystem in a rotor. Math- 
ematically, we use the Weyl pair (e*"^, V) of the rotor to 
define the Weyl pair {Z,X) of a qubit. Specifically, we 
choose the two unitary and hermitian operators Z and 
a: as 

Z = e"^ = (-1)^, 

A ^ + Z)V^ + V{1 + Z)). (11) 

They obey the relation of pU)) and hence they form the 
Weyl pair of a qubit. Note, furthermore, that the trace of 
Z and the trace of X vanish, as in the 2-dimensional case. 
The action of the operator Z is simple to understand as 
illustrated in Table [T] It only considers the parity of 
the eigenstate The total Hilbert space is split into 
two eigenspaces, one with eigenvalue -1-1 corresponding 
to even angular momenta and one with eigenvalue — 1 
for the odd angular momenta. Since the operator Z is 
infinitely degenerate, at least another operator is neces- 
sary to obtain a Complete Set of Commuting Operators 
(CSCO) of the original rotor. Since two adjacent angu- 
lar momenta are already discriminated by Z, the second 
operator of the CSCO can view them as identical. This 
can be done with the operator Li = [i/2j, where \x\ 
denotes the floor function, that is, the largest integer 
that does not exceed x. This is shown in Table [TJ Here, 
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we attach the index 1 to discriminate between the origi- 
nal rotor and what wih be found to be a residual rotor. 
The two operators Z and L\ commute and one can write 
£ = p + 2 [£/2j where [£/2j is the eigenvalue of Li and 
p such that z = (—1)'' where z is the eigenvalue of Z. 
Thus Z and L\ form a CSCO of the original rotor and 
completely specify the state of the system. We now have 
two alternative bases for the Hilbert space of the original 
rotor, either the eigenbasis of the operator L or the 
eigenbasis {|p,£i)} of the CSCO {Z,L-{). It remains to 
identify these two independent degrees of freedom and 
to make sure that their manipulation does not affect the 
other degree of freedom. In other words, first, we have to 
identify the two Weyl pairs of the two degrees of freedom 
and, second, not only Z and L\ have to commute but 
also their corresponding Weyl partners. 

For the operator Z, the Weyl pair is that of a qubit. 
With respect to the operator Li, one can substitute it 
with the operator e*"-^^ as they share the same eigenbasis. 
This suggests us to consider the pair (e'"^i,VL), where 
V\ = V'^ . The commutation relation of these two unitary 
operators is 

and e^'^'^^^ = 1. They form the Weyl pair of a rotor. As 
expected, this rotor is 27r-periodic. Indeed one can in- 
troduce the eigenvalue 6i of Vi by Vi \p, 9i) = e^^^ |p, 
and verify that 

|O,0i> = ^(|e> + |e + 7r>) -|0,^^l + 2^), 

|l,^i> = ^(l^>-|^ + 'r>) =|l,^?i + 27r), (13) 

where 9 — 9i/2. However, as noted earlier for the Fourier 
series of a 27r-periodic function, the ket 16*) can be decom- 
posed as 

\9) = %9,)+e-^'^'^\l,9^), (14) 

where the kets |O,0i) and |l,0i) are 7r-periodic in 9. 

Finally one can compute the four commutators 
[Z,e*"-^i], [X,e*"-^i], and [X,Vi] and ascertain 

that they vanish. Consequently the rotor has been de- 
composed into a qubit subsystem and another rotor. It 
is crucial to notice the tensor product structure of this 
mapping which, as mentioned above, can be cast as 
rotor = qubit ® rotor. In mathematical terms this trans- 
lates as H ~ H , where T-L denotes the Hilbert space 
of the rotor. Here the qubit is encoded as a subsystem 
rather than as a subspace. The advantage of encoding in 
a subspace or in a subsystem may appear when consid- 
ering a peculiar error model, a specific task to perform 
or, of course, practical considerations. 

Now we can write the logical qubit basis states as 

10,^-0)= "0^1^) ll'V-i)- I]"hK>,(15) 

I even I odd 
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TABLE 2: Encoding of three qubits in a rotor. Here three 
qubits are encoded by repeatedly factorizing the rotor degree 
of freedom. The residual rotor is obtained by concatenating 
2^ adjacent angular momentum eigenstates. 



where the coefficient a^i and an are complex numbers 
such that even l"o£p = odd Wu? ^ = ^ but other- 
wise unrestricted. The ket is labeled with the parity of 
the eigenvalue of the Z operator while ipo and ■01 spec- 
ify the state of the residual rotor, here assumed to be 
in a pure state. There is a freedom in the choice of the 
residual rotor states or equivalently in the coefficients au^ 
i = 0, 1. This is a consequence of the infinitely degener- 
ate spectrum of the operator Z. We will use this freedom 
to our advantage. More specifically we will exploit the 
degeneracy in two ways: first to encode many qubits, sec- 
ond, to protect these qubits against errors in the angular 
position and momentum. The encoding of many qubits 
is investigated in the remainder of this section while their 
protection against errors in realistic realizations will be 
thoroughly studied in Sec. IIVI 



B. Encoding many qubits in a rotor 

A consequence of the tensor product structure is to 
allow the encoding of many more logical qubits. Indeed 
nothing prevents us from exploiting the same encoding 
on the residual rotor to obtain two qubits, then three and 
even more. It should be clear that the above encoding 
allows in principle to encode as many qubits as desired. 
The limitation will of course come from the physical im- 
plementation as the more qubits we encode, the larger 
the accessible angular momenta must be. 

We shall now pay due attention to the Weyl pairs of 
different qubits. First we add an index j, j = 1, • • • , A^, 
to the operators Z and X to specify which qubit is con- 
sidered. The encoding of three qubits is exemplified in 
Table [51 This new mapping can be compared to the 
standard binary mapping recalled in Table [3] For the 
binary mapping, the encoding of the first qubit is anti- 
symmetric with respect to £ = while the encoding of all 
the remaining qubits is symmetric with respect io i — Q. 
This is in strong contrast with the new proposed map- 
ping which is invariant under translation of 2" angular 
momenta, where n is the number of qubits. For exam- 
ple, the encoding presented in Table [5] in invariant under 
translation of 8 angular momenta. By construction we 
can write the respective Weyl operators Zj and Xj for 
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TABLE 3: Binary encoding of three qubits. This encoding 
differs from the encoding presented in Table [2] 



the jth qubit as 

Xj = \{{t + Z,)V^^'-' +V^'-\t + Z,)). (16) 

The total Hilbert space is here again divided into two 
eigenspaces, corresponding to even and odd values of 
I £/2 ^~^ I , respectively. Evidently, the single qubit in 
Sec. IIII Al corresponds to the first qubit j = 1. Further- 
more, the Weyl pair of the residual rotor after encoding 
iV qubits is given by (e*"L^/^"J , V'^"). To illustrate the 
mapping presented above, let us explicitly write down 
the computational basis of the two-qubit subspace: 



|00,Voo) 
|01,Voi> 

110, ^10) 

111, ^11) 



00 

E 



= J2 "01^1 + 4^) 

00 

- aioe\2 + M) 

00 



(17) 



1, for 



with the proper normalization X]£=- 
j,k — 0,1. As already mentioned, the freedom in the 
residual rotor states will be used to ensure protection 
against shifts in angular position and momentum. 

We finally complete the set of single qubit operations, 
which can be constructed from Zj and Xj , with an entan- 
gling two-qubit gate to obtain a universal set of quantum 
gates. We choose the phase-shift gate R as the entangling 
two-qubit operation. This unitary transformation can be 
simply expressed in the {|-^)} basis as 



R 



jk 



E 

i=-oc 

Rkj , 



{-lY^]\i){£\ 
(18) 



where the indexes j and k denote the two qubits 
concerned by the gate so that £j — \£/2^~^\ and 

4= 



IV. PROTECTION OF MANY QUBITS 

The qubits encoded in this manner are vulnerable to 
errors. In the case of a single qubit, the logical qubits 
|0) and |1) are encoded in even angular momenta and 
odd angular momenta, respectively. Therefore, a state 
|0) affected by a unit shift £ £ + 1 will not be distin- 
guishable from an unaffected state 1 1) , as can be seen for 
Eq. ([15]). Distinguishability between |0) and |1) is then 
lost and it is thus necessary to provide an encoding which 
protects the logical qubits against small errors. 

In a very enlightening presentation Gottesman et 
al. use the stabilizer formalism to protect a qudit (sub- 
space) embedded in a infinite-dimensional space. Here 
we exploit the same technique to protect not only one 
but many qubits (subsystems). Actually we will show 
in Sec. I VIII that our results also apply to qudits. The 
technique used in ^5] slightly differs from the common 
use of the stabilizer formalism where k logical qubits are 
encoded in K > k physical qubits [3, ■ But the basic 
ideas are the same. 

In a nutshell the stabilizer formalism works as follows. 
We consider a set of orthogonal states that we want to 
protect against errors and call this subspace the code 
space. Within this code space, each state or code word 
is identified by the eigenvalue of a suitable operator. 
Errors are represented by unitary transformations that 
move the code space to another subspace. The differ- 
ent subspaces corresponding to different errors have to 
be distinguishable to allow correction. Therefore, they 
must be orthogonal. A simple solution to impose this or- 
thogonality is to associate each subspace with a different 
eigenvalue of a suitably chosen set of unitary transfor- 
mations. Eigenspaces are orthogonal and can therefore 
be perfectly distinguished. This chosen set of unitary 
transformations is the so-called stabilizer. By conven- 
tion we identify the code space as the -|-1 eigenspace of 
the stabilizer. The errors will then move the code words 
from the -1-1 eigenspace to another, therefore orthogonal, 
eigenspace of the stabilizer. In this ideal case the error 
is identified by reading out the eigenvalues of the stabi- 
lizer. Correction is then performed by applying the in- 
verse operator of the identified error. The -1-1 eigenspace 
is stabilized and contains our protected states. 

Here we are interested in two types of errors: continu- 
ous drifts for the angular position and discrete shifts for 
the angular momentum. They can be written in terms 
of the shift operators as 



and ELie)^^ 



(19) 



where e is real and e is an integer. To detect these two 
types of errors, two commuting unitary transformations 
are required. They also take the form of shift operators. 
A n-dimensional subspace of a rotor can be protected 
against shifts in angular momentum and angular position 
using the stabilizer 2] 



and Sl 



„2i7r-2-L 



(20) 
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where m is a free parameter related to the maximum 
amount of correctable errors. Sg will be used to read 
off the amplitude e of the error in 9, while Sl will be 
used to read off the amplitude e of the error in L. The 
above stabilizer can correct drifts in angular position up 
to |e| < tt/to and shifts in momentum up to |e| < m/(2n). 
For simplicity we now define A9 > and AL > as the 
maximal correctable error in 9 and L. For a fixed n, 
an increase in m will lead to an increase of AL but a 
decrease of A9. A trade-off follows, depending on which 
type of error is more critical for a given implementation 
or application. The stabilized space is composed of n 
code words or protected states 



\k[n], 1, 1) = y/m 



E 



\m{k/n + e)), 



(21) 



where k is an integer. Note that the factor -y/m is for con- 
venience. We use the label k[n] to emphasize that only 
k modulo n is relevant as can be seen from the above 
definition. The label 1,1 refers to the stabilizer's eigen- 
values. Let us also note that m/n has to be an integer to 
ensure the existence of a n-dimensional stabilized space. 
Otherwise only the state |0[n], 1, 1) = V^J2eL-rx> 1"^^) 
among the n states \k[n], 1, 1) is well-defined. We denote 
this integer by r. In the angular momentum basis the 
code words are equally weighted superpositions of shifted 
states. The code words take the form of a finite sum in 
the angular position basis 



\k[n],l,l) 



^ 'in — ± 

^ g-2j7rfej7n 
. /Til ^—^ 



3=0 



-J), (22) 

TO ' 



where the phases are powers of the nth root of unity. 

Since we want to protect N qubits, we have to identify 
the n-dimensional protected subsystem with the space 
spanned by the N qubits. Technically, not only the di- 
mensions have to match, i.e., n — 2^ , but also the qubit 
operators Zj and Xj have to commute with the other 
qubits' operators {Z/. and Xj., k ^ j) as well as with the 
stabilizer. A suitable encoding is provided by the follow- 
ing pairs 

z, = (-i)Lv(2-V)j^ 



+ Zj)Vt2^"^'' + + Z,)). (23) 



The CSCO of the original rotor now is 
{Zi,--- , Zn, Sl, Se). To prove it, let us consider 
an alternative writing of the basis {|^)}. 

The eigenvalues of Zj are of the form (— 1)^^ where 
Pj = 0, 1. The eigenvalues of Sl are of the form e^*'^''/'" 
where g = 0, • • • ,r — 1 while the eigenvalues of \_L/m\ 
are [i/mj. The identity 



N 



2^' V + g + [£/m\ to. 



(24) 
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TABLE 4: CSCO for two qubits and m = 6. Here the quan- 
tum numbers pi, p2, q, and [i/m] of the operators Zi, Z2, Sl, 
and lL/m\ uniquely specify the angular momentum eigen- 
state \£). 



therefore, provides a relabeling of the kets \£) in terms of 
the quantum numbers pj, q, and \£/m\. It follows that 
the operators (Zi,-- - , Z^, Sl, Se) commute and form 
a CSCO of the original rotor. Since V"^ and lL/m\ 
span the same subspace, the set (Zi, • • • , Zjv, Sl, Se) is 
a CSCO of the rotor too. 

An example of such CSCO is illustrated in Table [4] 
for two qubits and m — 6. For an easier reading of this 
table, the eigenvalues of the operators Zj and Sl are 
labeled by the powers pj and q instead of the eigenvalues 
themselves. For instance the eigenvalues of Sl for two 
qubits and m — 6 are 1^ e^*'^/'^, e"^''^/^ but we chose to 
label them as 0, 1, 2. 

One must not fail to note that (e*"Li/™J ^ y™) consti- 
tutes the Weyl pair of a rotor so that the last degree of 
freedom corresponds to a residual rotor. Moreover Xj 
commute with e^^L^-Z^J and V"^, with the other qubits' 
Weyl pairs as well as with Sl , so that we can finally write 



rotor — qubits (error in £) (g) rotor. 



(25) 



As a result we can still encode more qubits from the resid- 
ual rotor and the situation is similar to that of Sec. lIIII In 
principle we could also look at the Weyl pair associated 
with Sl and name the corresponding degree of freedom 
however this is irrelevant to the present study. 

The encoding of Eq. (j23p allows error correction in an- 
gular momentum up to AL < r/2. Thus we can choose 
to define r from the maximal number of correctable er- 
rors, that is, r = 2AL + 1. An even r is also conceivable 
however the resulting encoding would be less compact 
and will therefore require to access greater angular mo- 
menta. The simple encoding provided in Sec. lIIII is just a 
special case where AL = (r = 1) and would correspond 
to the stabilizer element Sl — ^ {m = n). The general 
protected qubits are then given by 



|A:H,l,l) = y'2^(2AL+l) ^ \{2AL + l){k + 2^ £)) . 

i=-oc 

(26) 

Since the protected states are N protected qubits, it is 
useful to label them with the eigenvalues of the opera- 
tors Zj. In the following we drop the eigenvalues of the 
stabilizer to lighten the notation. Let us illustrate what 
we have learned with three examples. First we can con- 
sider a single qubit with no protection against shift in the 



6 



angular momentum, that is, A'^ = 1, AL = 0, A6' < 7r/2. 
The computational basis is then given by 

oo 

|0) = ^/2 ^ \2l), 

£^-oo 
oo 

|1) = V2 ^ |l + 2£). 



(27) 



Equivalently, these states can be written in the angular 
position basis as: 



|o) = -L(|o) + k)), 
|i> - ^(|0)- w). 



(28) 



For a single qubit with protection against a unit shift in 
angular momentum (AL = 1 and A0 < 7r/6), we obtain 
the two protected states 



oo 

|1) = V6 ^ |3 + 6£). 



(29) 



A more interesting instance is concerned with two qubits. 
If we ask for protection against unit shifts in angular 
momentum (AL ~ 1 and A9 < 7r/12), we end up with 

oo 

|00) = 2V3 ^ 1 12^) , 
e=-oo 

oo 

|10) =2^3^ |3 + 12£), 

l=-oo 
oo 

|01) = 2^3 ^ |6 + 12£) , 



111) = 2^3 |9+ 12£) 



(30) 



As already mentioned earlier, these ideal code words 
are unphysical. This can be seen for example in Eq. (j27p 
and Eq. ([25)) . These states take the form of an infinite 
sum of equally weighted angular momenta and are per- 
fectly localized in angular position. In the following sec- 
tion we will investigate physical approximations of the 
ideal protected qubits. 



V. PHYSICAL APPROXIMATIONS 

The code words, and in particular the protected qubits, 
are infinitely squeezed states in angular position. With 
implementation in mind we should investigate the behav- 
ior of realistic approximations of the protected qubits. 



These are the physical approximations. Clearly, finite 
squeezing will inevitably lead to additional errors. Some- 
times these errors are within the range of correctable er- 
rors, sometimes they are not. We calculate the probabil- 
ity of noncorrectable errors for physical approximations 
of the protected qubits. 

We consider four physical approximations. Two ap- 
proximations pattern a finite squeezing in position. The 
first one relies on a truncated Gaussian wave function 
ll| while the second is a power of cosine wave function. 
The remaining two approximations are better envisioned 
in the angular momentum basis. One approximation re- 
sorts on a Gaussian envelop for the infinite sum of angu- 
lar momenta whereas the last instance focuses on a finite 
sum which we name grating. We study analytically the 
truncated Gaussian approximation and provide numeri- 
cal plots for the three other approximations. 



A. Error probability 

We want to define the probability of noncorrectable 
errors. For simplicity we first restrict ourselves to a sin- 
gle qubit only protected against drifts in angular posi- 
tion (AL = and A9 = tt/2). But the generality of the 
arguments remains and the formulas derived below are 
immediately generalized to larger N and AL. In view of 
Eq. (PSI) . the errors are specified on the two basis states 
|0) and Itt). Moreover measurements are here performed 
modulo TT in order to avoid destroying the superposition. 
Since only an error smaller than A0 can be corrected 
while the approximated basis state spreads all over the 
range [0,27r], some noncorrectable errors happen. For 
the state |0), any value of the angular position measured 
in the range [— 7r/2,7r/2] can be corrected while values 
within [— TT, — 7r/2] and [7r/2,7r] cannot. For the state |7r), 
the situation is opposite. Any value of the angular po- 
sition measured within [— tt, — 7r/2] and [k/2,it] can be 
corrected while values in [— 7r/2,7r/2] cannot. This nat- 
urally leads to the classical probability of error for the 
approximated states |0) and jvr): 



Pe{0) 
Pe(7r) 



do\{e\Q)\' 

d9\{e\7:)\\ 



de\{9\o) 



(31) 



In these equations and in the succeeding ones, unless 
stated otherwise, the symbols |0) and |7r) refer to the 
approximation of the basis states. Taking into account 
the symmetry of the wave function, we end up with 



Pe{0) =Pe(7r) 



d9 \ {9\tt) 



(32) 



We can now define the probability of error for any state 
diagonal in the basis {|0) , |7r)}. For a given density ma- 
trix p = q;o|0)(0| + Q;7r|7r)(7r| with ao + a-^ = 1, the error 
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probability is 

Pe{p) = aoPe{0) + a^Pe{T^)- 

This immediately leads to 



Pe{p) = 2 de\{e\TT} 

Jo 



(33) 



(34) 



Any equally-weighted classical mixture of the logical 
states |0) = (|0) + \n))/V2 and |1) = (|0) - \n))/V2 is 
diagonal in {|0) , |7r)}. Indeed, the off-diagonal elements 
nicely vanish. 

However, we cannot apply the above definition to more 
general mixtures (ao and ai different from one half) 
where the off-diagonal terms |0)(7r| and |7r)(0| do not van- 
ish. Nevertheless since the error is defined in the {\0)} 
basis we can reasonably extend the definition of the error 
probability to states diagonal in the {\0)} basis. 

Any good approximate state should have a sharp peak 
and fast decreasing tails. Therefore, a well approximated 
state should have nearly vanishing cross terms {9 \ 0){tt\ 9) 
and {9 |7r) (0| 6) leading to an almost diagonal state in the 
{\9)} basis. Under this constraint, the error probability 
for any approximated states is still given by Eq. ([M)) . 

The generalization to any number of qubits and any 
AL is immediately given by 



Peip) 



d9 \{9\Tr) 



(35) 



In the general case, the measurement of the angular po- 
sition is modulo 2Tr/m to avoid destroying the superposi- 
tion of Eq. . We compute the probability of error for 
the so-called truncated Gaussian states and numerically 
evaluate it for the three other approximations. 



B. Truncated Gaussian states 

The wave function of a truncated Gaussian state cen- 
tered at 9o and with a degree of squeezing ^ is given by 



(36) 



limit of no squeezing the truncated Gaussian wave func- 
tion reduces to the unit constant. Moreover the over- 
lap (0|O)(7r|0) decreases as ^ increases. Together with 
Eq. (1351) these properties imply that the error probabil- 
ity vanishes in the limit of infinite squeezing while in the 
limit of no squeezing, the error probability tends to the 
probability of a pure guess, that is, 1 — 1/(2^(2AL -I- 1)). 

Using Eqs. ([55)) - (|55| . the probability of error for the 
truncated Gaussian reads 



P. 



Iqubit 



1 



erf«/2) 
erf(7rO 



(38) 



for a single qubit. More generally, the error probability 
for N qubits is 



P, 



N qubits 



erf(7re/2^(2AL-H 1)) 



erf(7r^) 



(39) 



When ^ is of the order of m = 2 , the error probabil- 
ity qubits close to 10^^ for arbitrary N . In the 
limit of large squeezing, i.e., ^ ^ m, qubits g^^les as 
TOe~*^^^^/(7r'^/^^). The error probability is plotted in 
terms of the degree of squeezing ^ in Fig. [l]for one, five, 
and twenty qubits and AL — 1. 



Pe 
1.0 k 



AL=1 



0.8 




Log2<J 



FIG. 1: Error probability for the truncated Gaussian in the 
case of 1, 5, and 20 encoded qubits and protection against a 
single shift in angular momentum. When the physical param- 
eter ^ is much larger than the dimension of the qubits sub- 
space 2^, the error probability decays faster than a Gaussian. 
When ^ and 2^ are of the same order, the error probability 
is close to 10~^. 



where 6*0 is in [0, 27r] and 9q — 'k<9<9q + 'k. The nor- 
malization is = erf(7r^) where 



erf(a;) = —j= 



dt : 



(37) 



denotes the error function. As required, in the limit of in- 
finite squeezing the non-normalized wave function tends 
to 27r times the delta function. Note that formally the 
truncated Gaussian function is not 27r-periodic unlike the 
braket (6'|0o) = 2'n5^'^'"^{9 — 9q) however we only consider 
this function in the range Qq — ■k<9<9q^'k. In the 



C. Other approximations 

Here we consider the three remaining approximations. 
We only give their definitions and plot their error prob- 
abilities in terms of the relevant physical parameter for 
one and five qubits and Ai = 1. 

First the power of cosine wave function is defined as 




(40) 
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with the suitable normalization C^. Herep^ qubits close 
to 1 when 7 and m are of the same order. Therefore, the 
power 7 should be much larger than m to keep the level of 
error reasonably low. The graphs of the error probability 
in terms of the power 7 for one and five qubits are shown 
in Fig. El 



AL=1 




LogjY 



FIG. 2: Error probability for the power of cosine in the case of 
1 and 5 encoded qubits and protection against a single shift in 
angular momentum. When the physical parameter 7 and the 
dimension of the qubits subspace 2^ are of the same order, 
the error probability is close to 1. A reasonable approximation 
therefore requires 7 ^ 2^. 



Second the error probability for the Gaussian envelop 

(41) 



where is the required normalization, is plotted in 
terms of the width a in Fig. [31 Similarly to the truncated 



Finally, the error probability for the grating 

^= y 



(42) 



is represented in terms of the number of slits Lm in Fig.|4l 
Note that the range of this plot starts at (i.e., Lm = 1) 



Pe 
l.Ot 



0.2 



AL=1 

• N=l 
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FIG. 4: Error probability for the grating in the case of 1 
and 5 encoded qubits and protection against a single shift in 
angular momentum. For a number of slits Lm of the order of 
the dimension of the qubits subspace 2^, the error probability 
is close to 10~^. Consequently, the number of slits must be 
much larger than 2^ to obtain a satisfactory approximation. 



and the spacing between two points is not constant since 
Lm has to be an integer. The error probability qubits 
tends to 10^^ for when Lm and m are of the same order. 
Consequently, the number of slits must be much larger 
than m to obtain a low probability of error. 



VI. IMPLEMENTATION 



Pe AL=1 




N=l 

N=5 
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FIG. 3: Error probability for the Gaussian envelop in the case 
of 1 and 5 encoded qubits and protection against a single shift 
in angular momentum. As for the truncated Gaussian, the 
error probability is close to 10~^ when the physical parameter 
a and the dimension of the qubits subspace 2^ axe of the same 
order. 

Gaussian case, the error probability for the Gaussian en- 
velop is close to 10~^ for a and m of the same order. 



We discuss a possible implementation of a quantum 
rotor in quantum optics and briefly mention two other 
realizations using atom optics and molecular physics. 

In quantum optics, the orbital angular momentum of 
light can be used as a quantum rotor. It has been 
shown that any light beam with an amplitude distribu- 
tion of the form u{r,9,z) = uo{r, z)e^^^ , where £ is an 
integer, carries a quantized orbital angular momentum 
around the beam axis z [l^]. Physically the distribution 
u{r, 9, z) = uo{r, z)e*^^ can be obtained using a Laguerre- 
Gaussian (LG) light beam. There, each photon carries an 
angular momentum £h. In this context, the angular mo- 
mentum is also called LG mode. With today's technol- 
ogy LG modes can be easily produced and manipulated 
[13I. [1^. For the proposed encoding, an infinite super- 
position of LG modes is required. Such a superposition 
can be experimentally realized, e.g., using a single LG 
mode state which propagation axis is suitably displaced 
p^ . Obviously, in any experimental scheme only a su- 
perposition of a finite number of modes is achieved and 
approximations, like the grating or the Gaussian envelop 



presented in the Sec. |Vl have to be taken into account. 
Any operation on the encoded qubits can then be reahzed 
with holograms, lenses, and linear optical elements and 
the measurements are carried out with well-established 
interferometric technologies. 

Let us first consider the encoding of a single qubit 
with AL — 0. In this simplest case, the states of the 
computational basis are superpositions of either even or 
odd modes. To address the even andj)dd modes we use 
the sorting technique introduced in [13] where a Mach- 
Zehnder interferometer with a Dove prism inserted into 
each arm is used, as shown in Fig. [S^. By fixing the 
relative angle between the two Dove prisms to 7r/2 and 
correctly adjusting the path length of the interferometer, 
one can ensure that even modes appear in one port while 
odd modes appear in the other port of the interferome- 
ter. The qubit operation, e.g., Z is realized by applying 
a TT-phase shift to the odd modes, while the qubit oper- 
ation X is realized by passing the even modes through a 
hologram which increases all the momenta by one unit, 
and passing the odd modes through a hologram which 
decreases all the momenta by one unit. Finally, after 
applying the desired operations, the two arms can be re- 
combined. 

To encode two qubits with AL — 0, we sort the beam 
into four ports, as illustrated in Fig. [5}d. The first port 
corresponds to orbital angular momenta modulo 4, the 
second port corresponds to 1 modulo 4, and similarly 
for the third and fourth ports. This sorting is done by 
cascading Mach-Zehnder interferometers (with a Dove 
prism inserted into each arm) and holograms. The first 
interferometer sorts between even and odd angular mo- 
menta. The sorted mode are then passed through a sec- 
ond stage where they are sorted further. In this second 
stage, the relative angle between the two Dove prisms is 
7r/4. Even modes directly go through the interferome- 
ter, while odd modes go through an interferometer sand- 
wiched between a unit-decreasing hologram and two unit- 
increasing holograms (one for each output arm). Once 
the modes are sorted according to their modes modulo 
4, desired operations on the first and the second qubits 
can be simultaneously performed. Moreover, we can ap- 
ply an entangling operations, such as a controlled-not 
gate. Further studies would be required to include the 
measurement of the stabilizer. Let us emphasize the pe- 
culiarity of the present approach with respect to two 
well known alternatives that might seem at first sight 
equivalent. To encode, say, two qubits in optical modes, 
one could use four orbital angular momenta or even 
four spatial modes. Mathematically, the three encodings 
are very different. Indeed, as already mentioned, the 
present encoding corresponds to a two-qubit subsystem, 
that is, rotor = qubit (g) qubit (g) rotor, while the encod- 
ing using four angular momenta is a two-qubit subspace, 
i.e., rotor = (qubit (8) qubit) © rest, and the four spatial 
modes form a qudit system, i.e., 4 — qubit qubit. 

We now mention two other implementations. An al- 
ternative to photons is the use of one or many cold neu- 




(b) Possible implementation 



(=1 mod 4 

A^=+l holograms 



- 3 mod 4 



FIG. 5: Possible experimental setup, (a) A Mach-Zehnder 
interferometer with two Dove prisms at each arm is used to 
sort angular momentum states. By fixing the relative angle 
between the prisms to a = 7r/2 one can sort even and odd 
angular momenta, as described in (b) To encode two 

qubits with AL = 0, we sort the beam into four ports where 
each port corresponds to orbital angular momenta modulo 4. 
The sorting is done with Mach-Zehnder interferometers with 
a Dove prism inserted into each arm and holograms. Once 
the modes are sorted, desired operations on the first and the 
second encoded qubits can be simultaneously performed. 



tral atoms in a circular magneto-optical trap [H, [T7| . 
A magnetic field and several laser are used to create a 
ring-shaped potential where the atoms sit. Then the 
atoms are manipulated with laser pulses. Another op- 
tion is concerned with ultra cold molecules, typically 
carbon-based molecules with a quantized rotational mo- 
tion [l^, [l^. Promising molecules are polar diatomic 
molecules trapped in an optical lattice or a solid matrix 

MM- 

The three directions cited above are three possible op- 
tions for implementing a rotor operating in the quantum 
regime. They would deserve additional investigations. 



VII. ENCODING QUDITS IN A ROTOR 

In Sec, mil and IIVI we have presented a systematic way 
to encode and protect N qubits in a rotor. This con- 
struction can immediately be generalized to qudits. In 
this section we briefly summarize the relevant results. 

We can define the dth root of unity uj = e^"/'', where 
d stands for the dimension of the qudit. The Weyl pair 
{Z^j'^\ X^^^), j = 1, • • • ,N, where N specifies the number 
of qudits, takes the form Q 



X 



id) 



V 
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where P„(d) , is the projector onto the eigenspace 

of Z^^^ with eigenvalue +1. One can check that 
we stih have the required tensor product structure 
rotor = qudits ® (error va. rotor between the en- 
coded qudits, the errors in the angular momentum and 
the residual rotor. More specifically the Weyl pair of the 
residual rotor is (e»Li/(d"(2AL+i))j ^ ^d~(2AL+i)) 
the commutators between the Weyl pairs of the qubits, 
Sh and the Weyl pair of the residual rotor vanish. 

The stabilizer is Sl = e^*'^"^/™ and Se = V"™, where 
now n = and m = {2AL + 1). The protected qu- 
dits are still given by Eq. (gTI) and Eq. 

VIII. CONCLUSION 

We have presented a scheme to encode many genuine 
qubit subsystems in a rotor. We have shown how to 
manipulate and entangle them. We have also consid- 
ered a quantum error-correcting code to protect the many 
qubits against small errors in angular position and mo- 
mentum. The whole scheme is actually generalizable to 
qudits of any finite dimension. Furthermore we have 
considered physical approximations of the ideal scheme 
and their consecutive error probabilities. We have then 
turned to implementations and proposed several direc- 
tions. 

Theoretically, the advantage of encoding many qubits 
in a single degree of freedom is to allow their individ- 
ual manipulation with a single unitary transformation. 
Thus complex entangled states like cluster states can in 
principle be realized with a single manipulation of the 
rotor. In this article, good care has been taken to en- 
code genuine qubits subsystems in a rotor. A benefit is 
the freedom to use the remaining rotor to extract new 
qubits. Another advantage, inherent to the tensor prod- 
uct structure, is the possibility to correct errors. However 



the distinction between genuine qubits and pseudo qubits 
may become blurred after implementation. Nonetheless 
one might wonder whether genuine qubits and pseudo 
qubits exhibit a difference not only in terms of error 
model but also in operational terms. In other words, 
is there any task that can be performed with genuine 
qubits but that cannot be preformed equally well, or not 
at all, with pseudo qubits? 

The challenges finally lie in the implementation. The 
scalability of the proposed scheme is directly related to 
the accuracy the rotor can be manipulated with. There- 
fore, with today's technology, it seems reasonable to 
reach a dozen of encoded qubits. This would already 
allow the simulation of small quantum systems. 

Let us conclude with a few comments on the orbital an- 
gular momentum of light. Unlike earlier proposals, the 
full continuous variable character of this degree of free- 
dom is used, resulting in many genuine qubit subsystems. 
Moreover, qubits encoded in light are interesting not only 
because they are manipulated solely with linear optical 
elements but also because they are flying qubits. Thus 
one might conceive a computation taking place during 
the travel time between two remote places. 
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